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ABSTRACT 

We  rejuvenate  the  Chiellini  integrability  method  to  generate  some  new  nonlinear  integrable  equations  and  find 
their  analytic  solutions.  A  generalised  class  of  Van  der  Pol  equation  with  a  new  form  of  potential,  a  damped  Ermakov- 
Painleve  II  system  and  a  generalised  Milne-Pinney  equation  is  analysed. 
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1.  INTRODUCTION 

Ermakov  [1]  used  the  Milne  Pinney  equation  [2],  [3]  while  investigating  a  first  integral  for  the  corresponding 
time  dependent  harmonic  oscillator.  Since  then,  this  nonlinear  equation  has  gained  intensive  attention  [4],  [5]  of 
physicists  and  engineers  due  to  its  widespread  application  in  many  physical  problems  such  as  propagation  of  laser 
beams  in  nonlinear  media,  plasma  dynamics  etc.  It  is  well  known  that  the  general  solution  for  the  Milne-Ermakov- 
Pinney  equation 

y  +  (02{t)y  =  -^~,  (1.1) 

y 

where  y  denotes  double  differentiation  of  y  with  respect  to  time  t ,  (O—  0)(t)  is  a  time-dependent 

2  2  1/2 

frequency  function  and  K  is  a  numerical  constant  can  be  written  as  y  —  (Ax\  +  2Bxp2  )  ",  where  A.  B  and 

2 

C  are  constants  such  that  AC  —  B~  =  K  and  Xj  and  x2  are  two  independent  solutions  for  the  time-dependent 

2 

harmonic  oscillator  equation  X  +  C0~(t)x  =  0  . 

Equation  (1)  does  not  include  any  mechanism  of  damping.  Hence  it  is  natural  to  add  a  term  linear  in  the 
velocity,  yielding  the  damped  Milne-Ermakov-Pinney  equation 

y  +  /dy  +  (02(t)y--^r,  (1.2) 

y 


TRANS 

STELLAR 

•Journal  Publications  •  Research  Consultancy 


www.tjprc.ors 


editor@tjprc.  org 


2 


Supriya  Mukherjee,  A.  Ghose  Choudhury  &  Parthaguha 


where  jU  >  0  is  a  constant  positive  parameter.  Equation  (1.2)  can  be  transformed  into  generalized  Emden-Fowler 
equation  of  index  —  3  which  satisfies  integrability. 

In  recent  times  a  hybrid  Ermakov-Painleve  II  system  was  derived  by  Rogers  [6]  in  a  pioneering  work  as  a 
reduction  of  a  coupled  N  + 1  -dimensional  Manakov-type  NLS  system.  He  showed  that  the  Ermakov  invariants  admitted 
by  the  hybrid  system  were  key  to  its  systematic  reduction  in  terms  of  a  single  component  Ermakov-Painleve  II  equation 
which,  in  turn,  may  be  linked  to  the  integrable  Painleve  II  equation. 

The  force  free  Van  der  Pol  equation  and  the  Duffing  equation  are  characterised  by  the  equation 

x  =  g(x,x,t)  +  r(t)  with  g(x,x,t)  and  y(t)  intheform 

0  dV 

g(x,x,t)  =  (a-j3x  )x-——,  y(t)  =  0,  (1.3) 

dx 


1  2  2 

where  V  is  the  potential  function  approximated  by  a  finite  Taylor’s  expansion  in  series  V2(x)  —  and 


V4(x)  = —<VqX2  + —CHqX4  respectively,  OJ{)  and  Aq  are  non  zero  and  (X{]  is  a  real  number. 


The  Van  der  Pol  equation  introduced  by  Dutch  electrical  engineer  Balthazar  Van  der  Pol  in  1927,  has  gained 
considerable  importance  due  to  its  usefulness  in  modeling  various  physical  systems.  Since  its  introduction,  it  has  been  a 
prototype  for  oscillatory  systems  in  various  fields  of  science.  For  /U  =  0  the  equation  reduces  to  the  well  known  simple 

2 

harmonic  motion.  But  for  ±  0  the  term  //(I  —  x  )x  plays  an  important  role  in  describing  nonlinear  dissipative 
systems. 

The  Duffing  equation,  named  after  George  Duffing  in  1918  attracted  the  attention  of  many  researchers  due  to  its 
chaotic  nature  and  the  ability  to  replicate  similar  dynamics  in  our  natural  world.  This,  together  with  Van  der  Pol  equation 
has  become  one  of  the  most  common  examples  of  nonlinear  oscillations. 

Another  well  studied  but  still  challenging  equation  in  nonlinear  dynamics  is  the  combination  of  the  Duffing  and 
the  Van  der  Pol  equation.  The  nonlinearity  of  the  equation  makes  it  difficult  to  find  the  first  integrals  and  its  solutions.  Not 
much  progress  was  made  on  the  integrability  of  the  Duffing- Van  der  Pol  nonlinear  oscillator  equation  until  Chandrasekar 
et.al  [15]  have  established  complete  integrability  of  this  equation  and  have  derived  a  general  solution  for  a  specific  choice 
of  the  arbitrary  parameters. 

In  this  paper  we  have  introduced  a  generalised  Van  der  Pol  equation  with  a  new  form  of  potential  which  exhibits 
integrability.  An  analytic  solution  of  the  above  mentioned  system  is  also  derived  using  the  Chiellini  integrability  condition. 

The  application  of  the  Chiellini  integrability  condition  to  find  the  solutions  of  nonlinear  differential  equations  has 
been  recently  promoted  by  two  groups;  Harko,  Mak  and  their  coauthors  [16,  17,  18,  20]  and  Mancas  and  Rosu  [21,  22,  23]. 
It  must  be  worth  to  note  that  the  Chiellini  integrability  condition  appears  quite  naturally  for  Hamiltonization  of  the  Lienard 
equation  using  the  Jacobi  multiplier  technique.  Our  intention  in  this  letter  is  to  extend  the  scope  of  the  integrability 
condition  given  by  Chiellini  to  find  out  solutions  to  some  new  kind  of  nonlinear  systems.  In  particular,  we  obtain  the 
analytic  solutions  of  the  damped  Ermakov-Painleve  II  equation,  generalized  damped  Milne-Pinney  equation  and 
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generalised  Van  der  Pol  equation. 

2.  CHIELLINI  INTEGRABILITY  CONDITION  AND  ITS  APPLICATION 

The  first  order  Abel  differential  equation  [7]of  the  first  kind  plays  an  important  role  in  many  physical  and 
mathematical  problems.  The  connection  between  the  second-order  nonlinear  differential  equations  and  the  Abel  equation  is 
well  known  [8], [9]  and  the  solutions  to  such  differential  equations  can  often  be  obtained  via  the  solutions  of  the 
corresponding  Abel  differential  equations.  A  second  order  differential  equation  of  the  Lie'nard  type  [10]  given  by 

Y  +  g(Y)Y  +  h(Y)  =  0  (2.1) 

may  be  transformed  into  a  first-order  Abel  differential  equation  of  second  kind,  namely 

Z^~+g(Y)Z  +  h(Y)  =  0  (2.2) 

dY 

by  the  transformation  Y  =  z(Y (t))  which  in  turn  is  transformed  to  the  Abel  equation  of  first  kind 


—  =  g(Y)X2+h(Y)X3 


(2.3) 


Via  the  transformation  z.  =  —  ■  However,  the  criterion  of  integrability  of  such  equations  greatly  depends  on  the 

A 

expressions  of  g(Y)  and  h(Y).  An  important  observation  by  Chiellini  [11]  in  1931  states  that  a  first  kind  Abel  differential 
equation  (2.3)  is  exactly  integrable  if  the  functions  g(Y)  and  h(Y)  satisfies  the  condition 


d  f  h(Y)' 

dY{~g(Y), 


=  pg(X ) 


(2.4) 


For  some  constant  p.  This  integrability  condition  has  been  applied  in  1960s  by  Bandic'  who  wrote  a  couple  of 
mathematical  papers  [12], [13]  and  then  by  Borghero  and  Melis[14]  in  the  Szebehely’s  problem.  Recently,  this  integrability 
condition  as  gained  much  attention  by  Mak  and  Harko  [16], [17], [18]  in  obtaining  general  solutions  of  the  first-kind  Abel 
equations  from  a  particular  solution.  This  result  has  also  been  used  by  Yurov  and  Yurov  [19]  in  cosmology  and  again  by 
Harko  et  al  [20]  incase  of  particular  Lie'nard  equations. 

The  Chiellini  condition,  not  only  ensures  the  integrability  of  a  system,  but  it  also  helps  to  find  the  solution.  If  we 
h(Y) 

further  require  that  Z  —  ck  — yyy  then  its  substitution  in  equation  (2.2)  leads  to 


2  —l±Jl  —  4 ~p 

pck  +ck  + 1  =  0  =>  ck  = - 

2  p 

For  simplicity,  we  choose  ck=  1  which  gives  the  value  of  p  =  -2  in  equation  (2.4).  Thus,  from 
Y=  z(Y  (t))  we  have 
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Using  this  result  we  arrive  at  a  much  relevant  observation  that  equation  (2.1)  can  be  turned  to  the  non  dissipative 


equation 

Y  +  7/(7)  =  0,  H(Y)  =  21i(Y)  (2.6) 

where  the  function  h(Y )  is  scaled  up  by  a  factor  2  . 

This  result  allows  us  to  find  the  dissipation  function  in  (2.1)  without  actually  knowing  Y  .  Multiplying 
Y  +  2li(Y )  =  0  by  Y  and  integrating  we  have 


YY  +  2 Yh(Y)  =  0=>  F2  =  -A^li(Y)clY  +  < 


Where  c  is  an  integrating  constant.  Thus,  from  (2.5)  we  have 


$(r)=- 


%  J  h(Y)d(Y) 


Upon  further  integration  of  (2.7)we  have 


t~t 0  = 


%  J  h(Y)d{Y) 


where  (0  depends  on  initial  conditions. 


2.1.  Dissipative  Ermakov-  Painleve  II  Equation 

Combining  the  terms  of  both  Ermakov-Pinney  equation  and  the  Painleve  II  we  obtain  the  following  equation 


t  3  1  f  e 

F  +  +  =-— 3  r~- 

2  4 y3  v  2 


This  nonlinear  equation  is  known  as  the  (single  component)  Ermakov-  Painleve  II  equation  and  was  derived  by 
Rogers  etal  [24,6].  It  is  related  [25]  to  the  Painleve  II  equation 


z  =  2z  +tz+  y. 


(2.11) 


£  (  £ 

where  z  = - j\  7 - 2 yy 

2  y  V  2 


If  we  express  y  =  -\j\0\2  +\f/ 1“  then  the  canonical  single  component  Ermakov-  Painleve  II  equation  yields  a 

particular  Ermakov-Ray-Reid  system  (for  details,  see[6])and  admits  the  characteristic  invariant  which  may  be  exploited 
systematically  to  construct  the  solutions. 


Let  us  assume  the  equation  (2.10)  as 
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3  _  V 


y  +  Ay  +  ey  = 


(2.12) 


Where  A  =  —  and  77 
2 


4 


£ 

7 — 

V  2 


If  h(y)  =  Ay  +  £y3  -  ^ 

y 

then  equation  (2.12)  can  be  written  as 

yy  +  h(y)  =  0 


(2.13) 


(2.14) 


We  introduce  the  dissipative  Ermakov-  Painleve  II  equation  having  same  h(y)  as  in  the  non  dissipative  case  but 
with  an  additional  damping  term.  The  equation 


Y  +  g(Y)Y +  h(Y)  =  0 


(2.15) 


is  called  Chiellini  dissipative  Ermakov  Painleve  II  equation  because  the  damping  coefficient  g(Y  )  will  be 
obtained  from  the  Chiellini  integrability  condition.  Using  h(Y  )  from  (2.13)  in  (2.7),  we  have 


Y  =  -yJc-eY4-  2AY2  -  2 J]Y~2 

and  in  (2.8),  we  have 

!(f)= 

J-er6 -lAr*  +  cY2 -2tj 

Further  from  (2.16)  upon  integration  once  more,  we  have 


(2.16) 


(2.17) 


Y1  =\ 


{t  ) 

r  c2 


2+4 

\  1/2 


V 


lv16A2  Xj 


sin(2V2(7-f0  ))+■ 


4A 


:=  -1 


0 


(2.18) 


t0  depending  on  initial  conditions. 

2.2.  Generalized  Dissipative  Milne-Pinney  Equation 

At  first  we  embark  a  simple  equation  of  this  category  and  obtain  its  solution.  Let  us  consider  the  equation 

Y  +  g(Y)Y-4-  =  0 


(2.19) 


This  may  be  written  as  equation  (2.1)  with  h(Y)  = - —  . 
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Following  similar  arguments  we  have  Y  =  y  c  —  SY 


The  dissipative  term  g(Y)  = - ,  and  a  parametric  solution  to  the  equation  (2.20)  in  terms  of 

Y3  y  cY4  -  5 

Weierstrass  ^function  is  given  as 


2  p(co„)  16p(co„)  J  a(w+c-co0)  16  p((Oa)‘l 

{  P(a>0)  16  lP(ro„)  p(a>0)  Jl 


Y  =  yo+~ 


fiy o) 


4  £>(®+c)- 


/' (y0) 


where  J  is  an  integrating  constant  and  c  being  any  fixed  constant. 

Also  f(Y)  —  cY4  —  S,  y0  is  a  root  of  the  equation  /(F)  =  0  and  p  (ty)  =  p(69,g2,g^)  is  the 


Weierstrass  p 


function  attached  to  the  Weierstrass  Invariants.  Here  §2~  3Cf2  —  AcxfX^  and 


3  2  2  3 

g3  =  2flt1Cf2^3  ~(%2  ~  aoa\  !  (To  =  6' .  GCX-  cy0  ,  Cf2  =  (.'}'() ,  Cf3  =  cy0 .  <7(C0)  and  1]((0)  are  the  Weierstrass  sigma  and 

F 

/*  ^  ^ 

Weierstrass  zeta  functions  respectively,  p(<y0)  =  J — 04”  (for  a  choice  of  ry()  ). 

Generalizing  the  results  of  equation  (2.1),  we  consider  the  following  example  of  generalized  damped  Milne- 
Pinney  equation.  Let  us  consider  the  equation 

r  r  R 

Y  +  g(Y)Y+AY  =  ^  +  ^\+YJ$nY2n+'  (2.23) 

Y  Y  71=0 


For  R  =  0  we  have  F  =  ^2(S0  —  X)Y2  +c  —  4k2Y  l—2  k^Y 


(A-S0  )Y  -k2Y~2  ~kxY~3 
I^S0  -X)Y2  +  c- 4 k2Y~l  - 2 kxY~ 


The  dissipative  term  g(Y)  =  - 


Further  integration  of  (2.24)  yields  a  parametric  solution  in  terms  of  Weierstrass  ^function  given  as 


f  (y0)  1  <J(CO+C  +  OJ0)  _  „ 

t  =  y0a>+  '  -  -  log— - - — \+2(o)+c)C(o^)  +S 

4p(<y0)  <j(co+c-  co0) 


F  =  y0+- 


f'(y0) 


4  p(co+c)- 


/  (To) 
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Where  5  is  an  integrating  constant  and  c  being  any  fixed  constant. 

Also  /(7)  =  2(<Sq  -  A)Y  4  +cY 2  —  4k2Y  —  2k]  ,  y0  is  a  root  of  the  equation  f{Y)  =  0  and 
kp(CO)-kp(CO,g2,g 3)  is  the  Weierstrass p  -  function.  Here  g2  =  3oc2  —4(X\(X2  and  g3  =  2a{a2a3  -  0t2  -  a0a2', 


a0  =  2(S0  -  A) 

CCy  =  2(<5q  —  A)  Jq 


a2  -  2(Sq  -A)y 0  +  — 

6 

a3=2(S0-A)y30+3^-k2 

6 


■ft  \ 

jp(co0)  =  - —  (for  a  choice  of  0)Q  )is  not  equal  to  any  of  the  roots  of  4y3  -  g3y  ~  £3  =0 


2.3.  Generalised  VanDer  Pol  Equation  with  New  Form  of  Potential 

Let  us  consider  the  nonlinear  differential  equation 

Y  +  [a  +  pY2  +  yY)Y  +  h(Y)  =  0 

This  equation  may  be  written  as  Y  +  g(Y)Y  +  h(Y)  =  0 

Where,  g(Y )  =  a  +  jiY 2  +  yY  and  h(Y )  is  a  general  function  of  Y  (t).  Integrating  (2.4)  with 
p  =  -2  we  have/7(T)  =  g(Y)  —  2j*  g(Y)dY  +  c . 


(2.27) 

(2.28) 


(2.29) 


Where  c  is  an  integrating  constant.  From  (2.29)  for 


g(Y)  =  CC  +  pY2  +  ?Y  we  have 


h(Y)  =  J32Y5  + 1 fyY4  - (r  +  |)y3  +  (cjB - 3 ay)  +  (cy-2a2)Y  +  ca 


(2.30) 


dY  1  2  -x 

And  from  (2.5)  we  have  - =  c  —  2aY  —  y¥~ - 7 

dt  3 


(2.31) 


Integrating  once  we  have 


2  P 


3(r!  +<s3)  +  — (r  +  <f)  +  — -  — 
2/3  4  p2p 


[In 


r-r-<t>  +  — 
2£ 


■Jr- 


+  BY  +  C 


3(r  +  4>) 


3(r2  +  3>3)  +  ^(r  +  <f)-^A  +  —  3(r!  +  <s:)+^(r  +  <t)-^A  +  — 

p  ip1  p  p  ip2  p 


(2.32) 


2Y  +  B 


where  5  =  r+q>  +  Z, 

P 


c  =  r2  +<i>2  +^(r+o)+- 

P  P 


and 
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2r3  =  | 

24>3  = 

3.  CONCLUSIONS 

In  this  letter,  we  have  generated  some  second  order  nonlinear  damped  equations  using  the  Chiellini  condition  of 
integrability.  We  have  also  shown  that  they  are  exactly  integrable  and  have  computed  the  solutions.  If  the  coefficients  of  the 
second-order  nonlinear  equations  satisfy  some  specific  conditions  that  follow  from  the  Chiellini  integrabilty  condition,  then 
the  general  solution  of  such  equations  can  be  obtained  in  an  exact  parametric  form.  Moreover,  the  solution  to  the 
generalized  VanderPol  equation  are  obtained  for  arbitrary  values  of  a,  P  and  y. 
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